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Abstract 

Background: Trypanosoma cruzi, the agent of Chagas disease, is currently recognized as a complex of six lineages 
or Discrete Typing Units (DTU): Tcl-TcVI. Recent studies have identified a divergent group within Tel - TcI DO m- 
TcIdom. is associated with a significant proportion of human Tel infections in South America, largely absent from 
local wild mammals and vectors, yet closely related to sylvatic strains in North/Central America. Our aim was to 
examine hypotheses describing the origin of the TcI DO m genotype. We propose two possible scenarios: an 
emergence of TcI DO m in northern South America as a sister group of North American strain progenitors and 
dispersal among domestic transmission cycles, or an origin in North America, prior to dispersal back into South 
American domestic cycles. To provide further insight we undertook high resolution nuclear and mitochondrial 
genotyping of multiple Central American strains (from areas of Mexico and Guatemala) and included them in an 
analysis with other published data. 

Findings: Mitochondrial sequence and nuclear microsatellite data revealed a dine in genetic diversity across 
isolates grouped into three populations: South America, North/Central America and TcI DO m- As such, greatest 
diversity was observed in South America (A r = 4.851, n = 0.00712) and lowest in Tcl DOM (A r = 1.813, n = 0.00071). 
Nuclear genetic clustering (genetic distance based) analyses suggest that TcI DO m is nested within the North/Central 
American clade. 

Conclusions: Declining genetic diversity across the populations, and corresponding hierarchical clustering suggest 
that emergence of this important human genotype most likely occurred in North/Central America before moving 
southwards. These data are consistent with early patterns of human dispersal into South America. 

Keywords: Trypanosoma cruzi, Maxicircle, Microsatellite, Chagas Disease, Phylogeography, Population genetics, Tel 



Findings 

Trypanosoma cruzi, the aetiological agent of Chagas dis- 
ease, infects 6-8 million people in Latin America, while 
some 25 million more are at risk of acquiring the disease 
[1]. Parasite transmission to mammal hosts, including 
humans, can occur through contact with the faeces of 
hematophagous triatomine bugs. However, non-vectorial 
routes are also recognized, including blood transfusion, 
organ transplantation, congenital transmission, and oral 
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transmission via ingestion of meals contaminated with 
infected triatomine feces [2,3]. 

T. cruzi (family Trypanosomatidae; Euglenozoa: 
Kinetoplastida) is most closely related to several widely 
dispersed species of bat trypanosomes [4] . Salivarian try- 
panosomes including medically important Trypanosoma 
brucei subspecies, represent a more divergent group [5]. 
The age of the split between the T cmzZ-containing and 
T. bmcei-contzlrimg trypanosome lineages is thought to 
have been concurrent with the separation of Africa 
and South America/ Antarctica/ Australasia 100MYA [6], 
implying that T cruzi and the other Schizotrypanum 



© 2012 Zumaya-Estrada et al.; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the 
Creative Commons Attribution License (http://creativecommons.Org/licenses/by/2.0), which permits unrestricted use, 
distribution, and reproduction in any medium, provided the original work is properly cited. 



Zumaya-Estrada et al. Parasites & Vectors 201 2, 5:226 Page 2 of 9 

http://www.parasitesandvectors.eom/content/5/1/226 



Table 1 Trypanosoma cruzi 1 samples included in this study 


Strain 
code 


Strain 


Host/vector 


Country 


State Latitude 


Longitude 


Date 


Population 


Reference 


PALDA4 


PALDA4 


Didelphis 
olbiventris 


Argentina 


Chaco 


-27.133 


-61.460 


2001 


SOUTH 


Messenger et al., [12] 


PALDA21 


PALDA21 


Didelphis 
albiventris 


Argentina 


Chaco 


-27.1 33 


-61.460 


2001 


SOUTH 


Messenger et al., [12] 


PALDA5 


PALDA5 


Didelphis 
albiventris 


Argentina 


Chaco 


-27.133 


-61.460 


2001 


SOUTH 


Messenger et al., [12] 


PALDAV2 


PALDAV2A3 


Triatomo 
infestans 


Argentina 


Chaco 


-27.133 


-61.460 


2001 


SOUTH 


Messenger et al., [12] 


PALDA20 


PALDA20 


Didelphis 
olbiventris 


Argentina 


Chaco 


-27.133 


-61.460 


2001 


SOUTH 


Messenger et al., [12] 


COTMA38 


COTMA38 


Akodon 
boliviensis 


Bolivia 


Cotopachi 


-17.430 


-66.270 


2004 


SOUTH 


Messenger et al., [12] 


P234 


P234 


Homo sapiens 


Bolivia 


Cochabamba 


-17.380 


-66.160 


1985 


SOUTH 


Messenger et al., [12] 


P238 


P238 


Homo sapiens 


Bolivia 


Cochabamba 


-17.380 


-66.160 


1985 


SOUTH 


Messenger et al., [12] 


P268 


P268 


Homo sapiens 


Bolivia 


Cochabamba 


1 7 3QH 

- I /.ooU 


-66.160 


1987 


SOUTH 


Messenger et al., [12] 


SJM22 


SJM22 cl1 


Didelphis 
marsupialis 


Bolivia 


Beni 


-14.810 


-64.600 


2004 


SOUTH 


Messenger et al., [12] 


SJM34 


SJM34 


Didelphis 
marsupialis 


Bolivia 


Beni 


-14.810 


-64.600 


2004 


SOUTH 


Messenger et al., [12] 


SJM37 


SJM37 


Didelphis 
marsupialis 


Bolivia 


Beni 


-14.810 


-64.600 


2004 


SOUTH 


Messenger et al., [12] 


SJM39 


SJM39 cl3 


Didelphis 
marsupialis 


Bolivia 


Beni 


-14.810 


-64.600 


2004 


SOUTH 


Messenger et al., [12] 


SJM41 


SJM41 


Philander 
opossum 


Bolivia 


Beni 


-14.810 


-64.600 


2004 


SOUTH 


Messenger et al., [12] 


SJMC12 


SJMC12 


Philander 
opossum 


Bolivia 


Beni 


-14.810 


-64.600 


2004 


SOUTH 


Messenger et al., [12] 


XE5167 


XE5167 cl1 


Didelphis 
marsupialis 


Brasil 


Para 


-1.710 


-48.880 


1999 


SOUTH 


Messenger et al., [12] 


IM4810 


IM4810 


Didelphis 
marsupialis 


Brasil 


Manaus 


-3.070 


-60.160 


2002 


SOUTH 


Messenger et al., [12] 


B2085 


B2085 


Didelphis 
marsupialis 


Brasil 


Belem 


-1.360 


-48.360 


1991 


SOUTH 


Messenger et al., [12] 


XE2929 


XE2929 


Didelphis 
marsupialis 


Brasil 


Para 


-5.830 


-48.030 


1988 


SOUTH 


Messenger et al., [12] 


AAA1cl5 


AAA1cl5 


Rhodnius 
prolixus 


Colombia 


Casanare 


4.150 


-71.200 


2010 


SOUTH 


Ramirez et al., Molecular 
Ecology In press 


AAA7cl2 


AAA7cl2 


Rhodnius 
prolixus 


Colombia 


Casanare 


5.100 


-71.600 


2010 


SOUTH 


Ramirez et al., Molecular 
Ecology In press 


AAB3cl3 


AAB3cl3 


Rhodnius 
prolixus 


Colombia 


Casanare 


4.150 


-71.200 


2010 


SOUTH 


Ramirez et al., Molecular 
Ecology In press 


AAC1 cl3 


AAC1 cl3 


Rhodnius 
prolixus 


Colombia 


Casanare 


5.100 


-71.600 


2010 


SOUTH 


Ramirez et al., Molecular 
Ecology In press 


AACf1cl4 


AACf1cl4 


Canis familiaris 


Colombia 


Casanare 


5.100 


-71.600 


2010 


SOUTH 


Ramirez et al., Molecular 
Ecology In press 


AAD6cl6 


AAD6cl6 


Rhodnius 
prolixus 


Colombia 


Casanare 


5.100 


-71.600 


2010 


SOUTH 


Ramirez et al., Molecular 
Ecology In press 


CACQcl7 


CACQcl7 


Homo sapiens 


Colombia 


Santander 


6.963 


-73.420 


2009 


Tel DOM 


Ramirez et al., Molecular 
Ecology In press 


CACQcl8 


CACQcl8 


Homo sapiens 


Colombia 


Santander 


6.644 


-73.654 


2009 


Tel DOM 


Ramirez et al., Molecular 
Ecology In press 


DYRcl16 


DYRcl16 


Homo sapiens 


Colombia 


Boyaca 


5.640 


-72.899 


2007 


Tel DOM 


Ramirez et al., Molecular 



Ecology In press 
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Table 1 Trypanosoma cruzi I samples included in this study (Continued) 



EBclll EBclll 

FECcHO FECcHO 

Td3cl11 Td3d11 

X-1084cl10 X-1084cl10 

X-236cl9 X-236cl9 

YAS1cl3 YAS1cl3 

38 38 



Homo sapiens Colombia Boyaca 

Homo sapiens Colombia Boyaca 

Colombia Boyaca 

Colombia Boyaca 

Colombia Boyaca 



Triatoma 
dimidiata 



46 

66 

67 

70 

71 

83 

95 

100 

113 

116 

154 



CRISTY 

MICH1 

NINOA 

PLI 

QROI 

TQI 



Rhodnius 
prolixus 

Rhodnius 
prolixus 

Alouatta spp Colombia Casanare 



46 

66 

67 

70 

71 

83 

95 

100 

113 

116 

154 



Triatoma 
dimidiata 

Triatoma 
dimidiata 

Triatoma 
dimidiata 

Triatoma 
dimidiata 

Triatoma 
dimidiata 

Triatoma 
dimidiata 

Triatoma 
dimidiata 

Triatoma 
dimidiata 

Triatoma 
dimidiata 

Triatoma 
dimidiata 

Triatoma 
dimidiata 

Triatoma 
dimidiata 



DAVISc1 1 DAVIS 9.90 cl 1 Triatoma 
dimidiata 



ANITA II ANITA 

CAM6 CAM6 



Triatoma 
dimidiata 



CRISTY 
MICH 
NINOA 
PL 

QRO 
TQ 



Mexico 
Mexico 
Homo sapiens Mexico 
Mexico 
Homo sapiens Mexico 



Triatoma 
dimidiata 



Triatoma 
dimidiata 



Guatemala Jutiapa 

Guatemala Santa Rosa 

Guatemala Jalapa 

Guatemala Jutiapa 

Guatemala Jutiapa 

Guatemala Jalapa 

Guatemala Chiquimula 

Guatemala Chiquimula 

Guatemala Santa Rosa 

Guatemala Chiquimula 

Guatemala Baja Verapaz 

Guatemala Alta Verapaz 

Honduras Tegucigalpa 
Campeche 
Campeche 



Dipetalogaster Mexico 
maxima 



San Luis 
Potosi 

Michoacan 
Oaxaca 



Baja 

California Sur 



Triatoma 
barberi 

Triatoma 
pallidipennis 



Mexico Queretaro 
Mexico Morelos 



5.130 
5.920 
6.270 
4.960 
4.960 
5.300 
14.287 
14.177 
14.633 
14.287 
14.287 
14.633 
14.768 
14.768 
14.177 
14.768 
15.079 
15.594 
14.080 
19.188 
19.188 
22.159 
19.567 
1 7.054 
26.044 
20.594 
18.953 



-73.119 
-73.500 
-71.200 
-73.630 
-73.630 
-72.400 
-89.844 
-90.303 



-89.844 
-89.989 
-89.458 
-89.458 
-90.303 
-89.458 
-90.413 
-90.149 
-87.200 
-90.300 
-90.300 
-100.990 
-101.707 
-96.714 
-111.666 
-100.393 
-99.223 



2007 TclDOM 

2001 TclDOM 

2000 TclDOM 

2010 SOUTH 

2010 SOUTH 

2010 SOUTH 

2000 NORTH- 
CENT 

2001 NORTH- 
CENT 

2001 NORTH- 
CENT 

2001 NORTH- 
CENT 

2001 NORTH- 
CENT 

2001 NORTH- 
CENT 

2002 NORTH- 
CENT 

2002 NORTH- 
CENT 

2002 NORTH- 
CENT 

2002 NORTH- 
CENT 

2002 NORTH- 
CENT 

2002 NORTH- 
CENT 

1983 NORTH- 
CENT 

2011 NORTH- 
CENT 

2011 NORTH- 
CENT 

2007 NORTH- 
CENT 

2011 NORTH- 
CENT 

1994 NORTH- 
CENT 

2001 NORTH- 
CENT 

1986 NORTH- 
CENT 

1991 NORTH- 
CENT 



Ramirez et al., Molecular 
Ecology In press 

Ramirez et al., Molecular 
Ecology In press 

Ramirez et al., Molecular 
Ecology In press 

Ramirez et al., Molecular 
Ecology In press 

Ramirez et al., Molecular 
Ecology In press 

Ramirez et al., Molecular 
Ecology In press 

This study 
This study 
This study 
This study 
This study 
This study 
This study 
This study 
This study 
This study 
This study 
This study 

Messenger et al., 2012 

This study 

This study 

This study 

This study 

This study 

This study 

This study 

This study 
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Table 1 Trypanosoma cruzi I samples included in this study (Continued) 



XAL1 


XAL 


Triotomo 
dimidiata 


Mexico 


Veracruz 


19.173 


-96.133 


2003 


NORTH- 
CENT 


This study 




9209802P 


9209802P cl1 


Didelphis 
marsupialis 


USA 


Georgia 


32.430 


-83.310 


1992 


NORTH- 
CENT 


Messenger et al., 


[12] 


9307 


930701 03P ell 


Didelphis 
marsupialis 


USA 


Georgia 


32.430 


-83.310 


1993 


NORTH- 
CENT 


Messenger et al., 


[12] 


ARMA 


USAARMA cl3 


Dasypus 
novemcinctus 


USA 


Lousiana 


30.500 


-91.000 


Unknown 


NORTH- 
CENT 


Messenger et al., 


[12] 


USA 


USAOPOSSUM 
c\2 


Didelphis 
marsupialis 


USA 


Lousiana 


30.500 


-91.000 


Unknown 


NORTH- 
CENT 


Messenger et al., 


[12] 


9354 


9354 


Homo sapiens 


Venezuela 


Sucre 


10.460 


-63.610 


1999 


Tel DOM 


Messenger et al., 


[12] 


1 1541 


1 1541 


Homo sapiens 


Venezuela 


Merida 


8.590 


-71.230 


2003 


Tel DOM 


Messenger et al., 


[12] 


11713 


11713 


Homo sapiens 


Venezuela 


Lara 


10.233 


-69.866 


2003 


Tel DOM 


Messenger et al., 


[12] 


11804 


11804 


Homo sapiens 


Venezuela 


Portuguesa 


9.084 


-69.103 


2003 


Tel DOM 


Messenger et al., 


[12] 


1 0462P2C3 


1 0462P2C3 


Homo sapiens 


Venezuela 


Miranda 


1 0.266 


-66.485 


Unknown 


Tel DOM 


This study 




10462P2C7 


10462P2C7 


Homo sapiens 


Venezuela 


Miranda 


10.080 


-66.449 


Unknown 


Tel DOM 


This study 




10968P1C1 


10968P1C1 


Homo sapiens 


Venezuela 


Sucre 


10.406 


-63.298 


Unknown 


Tel DOM 


This study 




ANT3P1C6 


ANT3P1C6 


Homo sapiens 
(oral) 


Venezuela 


DC 


10.500 


-66.951 


Unknown 


SOUTH 


This study 




M13 


M13 


Didelphis 
marsupialis 


Venezuela 


Barinas 


7.500 


-71.230 


2004 


SOUTH 


Messenger et al., 


[12] 


M16 


M16 cl4 


Didelphis 
marsupialis 


Venezuela 


Barinas 


7.500 


-71.230 


2004 


SOUTH 


Messenger et al., 


[12] 


M18 


M18 


Didelphis 
marsupialis 


Venezuela 


Barinas 


7.500 


-71.230 


2004 


SOUTH 


Messenger et al., 


[12] 


M7 


M7 


Didelphis 
marsupialis 


Venezuela 


Barinas 


7.500 


-71.230 


2004 


SOUTH 


Messenger et al., 


[12] 


92122 


92122102R 


Procyon lotor 


TcIV 


USA 


Georgia 






OUTG ROUPS 


Messenger et al., 


[12] 


CANIII 


CANII cl1 


Homo sapiens 


TcIV 


Brazil 


Belem 






OUTGROUPS 


Messenger et al., 


[12] 


CM17 


CM17 


Dasypus spp. 


Tel II 


Colombia 


Carimaga 






OUTGROUPS 


Messenger et al., 


[12] 


X1060 


X10610 cl5 


Homo sapiens 


TcIV 


Venezuela 


Guarico 






OUTGROUPS 


Messenger et al., 


[12] 


ERA 


ERA cl2 


Homo sapiens 


TcIV 


Venezuela 


Anzoategui 






OUTGROUPS 


Messenger et al., 


[12] 


10R26 


10R26 


Aotus spp. 


TcIV 


Bolivia 


Santa Cruz 






OUTGROUPS 


Messenger et al., 


[12] 


SAIRI3 


Saimiri3 ell 


Saimiri 
sciureus 


TcIV 


Venezuela 


Venezuela 






OUTGROUPS 


Messenger et al., 


[12] 



species evolved exclusively in South America. Others 
propose an alternative origin of T. cruzi from an an- 
cestral bat trypanosome potentially capable of long 
range dispersal [7]. Whilst the precise scenario for the 
arrival of ancestral Schizotrypanum lineages in South 
America is a matter for debate, the current continental 
distribution and genetic diversity of T. cruzi supports 
an origin within South America. Parasite transmission 
is maintained via hundreds of mammal and triatomine 
species in different biomes throughout South and 
Central America, as well as the southern states of the 
USA [8]. 

Biochemical and molecular markers support the exist- 
ence of six lineages or Discrete Typing Units (DTU): Tel, - 



TcVI agreed by international consensus ([9]. Each DTU 
can be loosely associated with a particular ecological and/ 
or geographical framework [10]. Tel is ubiquitous among 
arboreal sylvatic foci throughout the geographic distribu- 
tion of T. cruzi and is the major agent of human Chagas 
disease in northern South America. Several molecular 
tools now identify substantial genetic diversity within Tel 
[11-14]. Importantly these new approaches consistently 
reveal the presence of a genetically divergent and homoge- 
neous Tel group (henceforth TcI DO m - previously Tela/ 
VEA/dqm) associated with human infections from Vene- 
zuela to Northern Argentina, and largely absent from wild 
mammals and vectors sampled to date [14]. The origin of 
this clade is unclear, although recent work supports a 
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C/D 
<f) 
(D 

O 

DC 
o 

< 

1— 
CO 
0) 

o 



Soi/fA? America 



North / Central 
America 



Teh 



DOM 



0.008 
0.007 
0.006 
0.005 
0.004 
0.003 
0.002 
0.001 



Figure 1 Nucleotide diversity and allelic richness comparisons across North and South American. Trypanosoma cruzi I populations. Left 
hand data points (diamond) indicate allelic richness ± standard error over loci. Right hand data points (square) indicate nucleotide diversity (n) : 
standard error over pair-wise comparisons. 



sister group relationship with Tel circulating in North 
America (e.g. [12,13]). 

In this manuscript we have set out to evaluate the genetic 
diversity of Tel in North/Central America, undertaking a 
comparison with Tel diversity in South America, including 
TcI DO m- Our aim was to examine hypotheses describing 
the origin of the TcI DO m clade. We propose two possible 
scenarios: an emergence of TcI DO m in northern South 
America as a sister group of North American strains and 
dispersal among domestic transmission cycles, or an origin 
in North America, prior to dispersal back into South 
American domestic cycles, possibly anthropically. To pro- 
vide further insight into this question we undertook high 
resolution nuclear and mitochondrial genotyping of mul- 
tiple Central American strains (from areas of Mexico and 
Guatemala) and included them in an analysis with other 
published data [11-13]. 

A panel of 72 Tel isolates and clones was assembled for 
analysis (Table 1) [11-16]. Of these, existing sequences 
and microsatellite data were available for 46 isolates 
[11,12]. Isolates were classified into three populations: Tcl- 

NORTH-CENT> TcI S oUTH an d TcI DO M- TcI NO RTH-CENT 

includes samples from the USA, Mexico, Guatemala and 
Honduras; TcI S outh corresponds to South America 
(Argentina, Bolivia, Colombia, Venezuela and Brazil) and 
TcI DO m with exclusively domestic isolates from Colombia 
and Venezuela, already known to correspond to a geno- 
type with restricted genetic diversity: Tela, as previously 
described by Herrera et al, (2007) [17] and VEN Dom , as 
described by Llewellyn et al, (2009) [13]. Additional DTU 
isolates (TcIII-TcIV) were included as out-groups in the 
mitochondrial analysis. 

Isolates from Mexico and Guatemala were character- 
ized to DTU level via the amplification and sequencing 



of glucose-6-phosphate isomerase (GPI) as previously 
described by Lauthier et al, (2012) [18]. Subsequently, 
nine maxicircle gene fragments were amplified, 
sequenced and concatenated from the Mexican and 
Guatemalan strains according to Messenger et al, 2012 
(excluding ND4) [12]. Phylogenetic analysis was also 
conducted as in Messenger et al, 2012 [12]. Nineteen 
nuclear microsatellite loci previously described by 
Llewellyn et al, 2009 [13], were selected based on their 
level of Tel intra-lineage resolution. Microsatellite loci 
were amplified across 21 unpublished biological stocks 
from Mexico and Guatemala. Reaction conditions were 
as described previously [13]. Dendrograms based on 
multilocus allele profiles were constructed also accord- 
ing to Llewellyn et al, 2009 [13]. 

Maxicircle nucleotide diversity (tt) was calculated for 
TcI NO rth-cent> TcIsouth and TcI D om respectively in 
DNAsp v5 [19]. Nuclear allelic diversity was calculated 
for the same populations using allelic richness (A r ) in 
FSTAT [20]. The resulting values are shown in Figure 1. 

Nucleotide sequences per gene fragment are available 
from GenBank under the accession numbers MURF1 
(fragment a): JX431060 - JX431084; MURF1 (fragment 
b): JX431085 - JX431109; ND1: JX431110 - JX431134; 
NDS (fragment a): JX431135 - JX431159; NDS (fragment 
b): JX431160 - JX431184; 9S rRNA: JX431185 - 
JX431209; 12S rRNA: JX431210 - JX431234; COIL 
JX431235 - JX431259 and CYT b: JX431260 - JX431284. 

Across the 3,449 bp final concatenated alignment (in- 
cluding outgroups), a total of 374 variable sites were 
found. The mitochondrial phylogeny supported the pres- 
ence of significant diversity among the isolates examined 
(Figure 2). TcI DO m strains formed a monophyletic clade 
[60% ML BS/0.98 BPP]. The principal division in the 
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83 
DAVIScIl 
. 9209802p 
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ANITA II 
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PL I « 



NORTH- 
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M16 
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CENT 
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0.97 
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SJM41 
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— COTMA38 
I P238 
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USA 




Mexico 


0 Intra Tel introgression 
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Honduras 






Venezuela/ Human 






Colombia/ Human 
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Figure 2 Isolate grouping of 72 Trypanosoma cruzi I strains, as well as outgroups, based on nine concatenated maxicircle sequences. 

Bayesian consensus topology is displayed. Bayesian posterior probability analysis (BPP) was performed using MrBAYES v3.1. Five independent 
analyses were run using a random starting tree with three heated chains and one cold chain over 10 million generations with sampling every 10 
simulations (25% burn-in). Decimal values (second number) on nodes indicate Bayesian probabilities for clusters. First number indicates the 
Maximum-Likelihood (ML) % bootstrap support for clade topologies, which was estimated following the generation of 1000 pseudo-replicate 
datasets. Branch colours indicate isolate origin. Isolates that show clear incongruity between nuclear genotype and maxicircle genotype are 
marked. Outgroup branches were cropped for ease of visualization, full branch lengths are show inset top right. 



phylogeny was between TcI SO uth and TcI DO m/TcI- 
north-cent (98% ML BS/0.98 BPP). However, this div- 
ision is incomplete, such that a subset of South 
American strains is also grouped with TcI DO m and Tcl- 
north-cent- Thus, it is not possible to conclude that 
TcI DO m maxicircle sequences nest uniquely among those 
from TcI NO rth-cent strains. Conversely, a basal rela- 
tionship of the TcI NO rth-cent to TcI DO m is suggested at 
the level of nucleotide diversity by population (Figure 1), 
whereby TcI D om<TcI NO rth-cent<TcI S outh- Low stand- 
ard errors about the mean in all three populations, but 
especially in TcI DOM and TcI NO rth-cent> suggest that 
sample size had little impact on the accuracy of estima- 
tion between populations. 



Distance-based clustering using the microsatellite data- 
set indicated the presence of several well defined clades 
(Figure 3). Importantly in this case the monophyly of 
North-Central American isolates was corroborated, and 
TcI DO M clustered firmly within them (bootstrap 65%). By 
contrast, South American isolates fall into a divergent but 
diverse clade. Thus the nuclear data provide stronger sup- 
port for divergence of TcI DO m fr° m within TcI NO rth-cent 
than the maxicircle phylogeny. Sample size-corrected al- 
lelic richness estimates are consistent with hierarchical 
patterns of clustering based on pair-wise genetic distances. 
As with the maxicircle dataset, there is a pronounced cline 
in diversity across the populations studied - A r TcI DO m< 
A r TcI NO rth-cent< A r Tclsouth (Figure 1). 
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Figure 3 Isolate grouping of 72 Trypanosoma cruzi I strains based on nineteen nuclear micrsoatellite markers. Neig hour-joining 
clustering algorithm implemented. Bootstrap values are included on important nodes. The first figure indicates % bootstrap support over 10,000 
trees, the second the % stability over 1000 trees accounting for multi-allelic loci in the dataset. For further details see Llewellyn et al., 2009 [13]. 
Branch colours indicate isolate origin. The three principal populations TcI DO m TcI SO uth and TcI NO rth-cent are shown on both map and tree. Red 
circles correspond to isolates from TcI DO m- Isolates that show clear incongruity between nuclear genotype and maxicircle genotype are marked 
with reference to Figure 2. 



Tel dispersion into Central and North America 

Using a 100 MYA biogeographic calibration point [6], 
molecular clock analyses point to the origin of T. cruzi 
(sensu stricto) 5-1 MYA [21-23] and a most recent 
common ancestor for Tel at 1.3-0.2 MYA [22]. Reduced 
genetic diversity among North-Central American isolates 
by comparison to their southern counterparts is power- 
ful evidence in support of others who suggest that Tel 
originated in South America [13,24]. The emergence of 
Tel in the South occurred prior to either migration 
across the Isthmus of Panama alongside didelphid mar- 
supials during the Great American Interchange [25], or 
perhaps prior to northerly dispersal via volant mammals 
(e.g. bats). 

Origin of TcI DO m 

Recent findings indicate a close resemblance between 
TcI DO m isolates from the northern region of South 



America and parasite populations from Central and 
North America by the use of nuclear and mitochondrial 
markers [11-13]. Indeed SL-IR genotyping suggests a 
distribution for TcI DO m that now extends as far south as 
the Argentine Chaco, where multiple sequences have 
been identified from human and domestic vector 
sources [14]. Llewellyn et al, (2009) originally hypothe- 
sised that a distinct human/domestic clade could be 
maintained despite the presence of nearby infective syl- 
vatic strains due to the low parasite transmission effi- 
ciency by the vector [13]. In this case multiple feeds by 
domestic vector nymphs are required to infect indivi- 
duals, as such human - human transmission is far more 
common than reservoir host - human transmission. Ori- 
ginally this hypothesis was developed to explain the epi- 
demiology of Chagas disease in Venezuela. However, 
TcI DO m is clearly widespread and recent data propose a 
date for its emergence 23,000 ± 12,000 years ago [11]. 
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This corresponds to the earliest human colonisation of 
the Americas [26]. Thus it seems that TcI DO m may be as 
ancient as humans in South America. Crucially, our 
data, which show that TcI DO m is nested among North 
and Central American strains, suggest that this wide- 
spread domestic T. cruzi genotype may actually have 
made first contact with man in North- Central America. 

The expansion of limited diversity genotypes into do- 
mestic transmission cycles is a familiar story in T. cruzi. 
This phenomenon seems to have occurred almost simul- 
taneously with TcI DO m (< 60,000 YA) in the Southern 
Cone region but involving DTUs TcV and TcVI [22]. 
Nonetheless, static human population densities sufficient 
to support a sustained domestic cycle are presumably 
vital. For TcI DO m> patterns of genetic diversity suggest 
early colonizing Amerindians may have been responsible 
for its southerly migration and dispersal from North/ 
Central America. However, such early settler populations 
were probably small, dynamic, and inherently unsuitable 
to sustain transmission of such a genotype. Many ques- 
tions, therefore, remain unanswered regarding its emer- 
gence. Insight could perhaps be drawn from a better 
understanding of the current distribution and diversity 
of TcI DO m (including samples from the Southern Cone), 
patterns of vector population migrations, and even from 
analysis of ancient DNA (e.g. [27]). We hope this report 
serves to galvanize efforts towards this understanding, 
especially among researchers in Central and North 
America, where many of the answers lie. 
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